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Abstract 

We analyze soft and collinear gluon resummation effects at the N^LL level for Stan¬ 
dard Model Higgs boson production via gluon fusion gg ^ H and the neutral scalar 
and pseudoscalar Higgs bosons of the minimal supersymmetric extension at the 
N^LL and NNLL level, respectively. We introduce rehnements in the treatment of 
quark mass effects and subleading collinear gluon effects within the resummation. 
Soft and collinear gluon resummation effects amount to up to about 5% beyond the 
fixed-order results for scalar and pseudoscalar Higgs boson production. 
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1 Introduction 


The Standard Model (SM) of elementary particle physics has been established as a very 
successful theory that describes the properties and strong, weak and electromagnetic 
interactions of all known elementary particles [1]. The recent discovery of the Higgs boson 
with a mass of about 125 GeV [2] completed the particle content of the SM and established 
the Higgs mechanism for electroweak symmetry breaking [3]. Its couplings to SM gauge 
bosons, i.e. to ZZ^ hT+hT“, and fermion pairs (r leptons and bottom quarks) as well as 
the loop-induced couplings to gluon and photon pairs, have been measured with accuracies 
of 10-50%. All measured couplings are in agreement with the SM predictions within their 
uncertainties [1]. In addition there are very strong indications that the newly discovered 
boson carries zero spin and positive CP-parity [5]. Its discovery is of vital importance for 
the mathematical consistency of the SM and the success of the predictions for the precision 
electroweak observables which are in striking agreement with measurements at LEP and 
SLC [6]. Based on the present situation it is highly relevant to test the properties of the 
discovered particle in more detail. The measured inclusive production and decay rates 
are in agreement with the theoretical predictions within the corresponding uncertainties. 

At hadron colliders as the LHC neutral Higgs bosons of the Standard Model are 
copiously produced by the gluon fusion process gg —)■ iL, which is mediated by top and to 
a lesser extent bottom and charm quark loops. Due to the large top Yukawa coupling and 
the large gluon luminosities gluon fusion comprises the dominant Higgs boson production 
mechanism for the SM [7]. 

In the past the NLO QCD corrections to the top and bottom/charm quark loops have 
been calculated [811^101111] . They increase the cross sections by up to 90%. The full 
quark and Higgs mass dependencies for the quark loops have been included [9II10II11] . The 
scale dependence decreased from 0(100%) to ~ 20% at NLO. The NLO results indicated 
that the limit of heavy top quarks provides a reliable approximation of the full relative 
QCD corrections within ~ 5% for a Higgs mass of 125 GeV [9|T^ . In this limit the 
cross section factorizes into a part originating from an effective Lagrangian derived by 
integrating out the top quark supplemented by gluonic and light quark corrections within 
the effective low-energy theory M. Within this approach the NNLO QCD corrections 
have been calculated leading to a further increase of the cross section by 20-30% [mil5lll6| 
and a further reduction of the scale dependence to ~ 10%. Mass effects beyond the 
heavy top approximation have been studied at NNLO in a heavy top mass expansion 
HT]. These mass effects amount to less than a per cent in the relative QCD corrections 
below the tt-threshold. The NNLO results have been improved by a soft and collinear 
gluon resummation at the NNLL level [TSlIT^ that add another 5-10% beyond NNLO to 
the total Higgs production cross section. In addition the small quark mass effects have 
been included in the soft-gluon-resummed result at NLL [20] • Very recently the N^LO 
corrections have been computed in the limit of heavy top quarks first approximately 
deriving the leading terms in a threshold expansion [21] and then completely [22)123] . 
These increase the cross section by only a few per cent and reduce the scale dependence 
to a level of about 5%. They coincide with the soft-l-collinear gluon approximation at 


1 









N^LO at the level of a few per cent [2^251126] . This signalizes that soft gluon effects 
are less important at N^LL than at NNLL level thus underlining a proper perturbative 
reliability of the N^LO and N^LL results. These results have been completed by the 
calculation of the full electroweak corrections at NLO 123 and beyond [28] that provide 
a further increase of the cross section by about 5%. 

One of the most attractive extensions of the SM is the minimal supersymmetric SM 
(MSSM) that requires the introduction of two Higgs doublets and thus predicts the exis¬ 
tence of hve elementary Higgs bosons, two neutral CP-even h, H, one neutral CP-odd A 
and two charged ones [29]. A LO the Higgs sector is described by two independent in¬ 
put parameters that in the case of a real MSSM are usually chosen to be the pseudoscalar 
Higgs mass Ma and the parameter tg(3 dehned as the ratio of the two CP-even vacuum 
expectation values. Higher-order corrections to the MSSM Higgs masses and couplings 
turned out to be large with the dominant piece originating from contributions induced by 
the large top Yukawa coupling [30] . Including all relevant corrections up to the three-loop 
level the mass of the light scalar Higgs boson h is bound to be smaller than ~ 135 GeV 

[31] . The light scalar Higgs boson becomes SM-Higgs like close to its upper mass bound 
thus allowing for the possibility of the Higgs boson found at the LHC to be identihed 
with the light scalar h. Global hts to the MSSM leave a small region within the MSSM 
parameter space where the discovered Higgs particle could also be the heavy scalar H 

[32] but this region is disfavoured by LHC data [33] . The Higgs couplings to intermediate 
gauge bosons and fermions are modihed by additional factors shown in Table [1] that only 
depend on the mixing angles a and /5 where a denotes the mixing angle between the 
neutral CP-even Higgs states. For large values of tg(3 the down-type Yukawa couplings 
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Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to the SM couplings. 


are strongly enhanced, while the up-type Yukawa couplings are suppressed, unless the 
light (heavy) scalar Higgs mass ranges at its upper (lower) bound, where the couplings 
become SM-sized. 

The dominant Higgs boson production mechanisms for small and moderate values of 
tg/9 are the gluon fusion processes gg —)■ h,H,A that are mediated by top and bottom 
and in addition stop and sbottom loops in the case of the scalar Higgs bosons h, H. For 
large tg/3 the leading role is taken over by Higgs radiation off bottom quarks due to the 
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strongly enhanced bottom Yukawa couplings [3l] . In the past the NLO QCD corrections 
for pseudoscalar Higgs boson production have been derived in the limit of heavy top 
quarks [5S] and later including the full Higgs and quark mass dependencies pi55] . They 
are large ranging up to about 100%. In the limit of heavy top quarks the calculation of the 
SM Higgs boson at NNLO QCD has been extended to the pseudoscalar case, too [laaTj 
inducing a moderate increase of the total cross section by about 20-30% as in the SM 
Higgs case. However, the limit of heavy top quarks is only applicable for the pure top loop 
contributions so that for bottom-loop dominance for large values of tg/9 we are left with 
NLO accuracy. For the case of top-loop dominance, i.e. for small and moderate values 
of tg/3 within the MSSM, the soft-gluon resummation has been extended to the NNLL 
level for the pseudoscalar Higgs boson [38] resulting in corrections of 0(10%) beyond 
NNLO. For the light and heavy scalar MSSM Higgs bosons the NLO QCD corrections 
to the squark loops have been calculated in the approximate limit of heavy squarks [32] 
yielding NLO corrections close to 100% as for the quark loops. The next-to-leading order 
computation for the squark loops including the full mass dependencies has been performed 
later [TTMO] for the pure QCD corrections. The genuine supersymmetric QCD corrections 
have been derived in Ref. [TT] for large supersymmetric particle masses, i.e. the full Higgs 
mass dependence has not been taken into account. For gluino masses much larger than 
the stop and top masses the results of [H] develop a logarithmic singularity in the gluino 
mass that seems to contradict the Appelquist-Carazzone decoupling theorem [12| at hrst 
sight. This problem has been solved by a detailed renormalization group analysis with 
the corresponding decoupling techniques for heavy particles [33|. The full supersymmetric 
QCD corrections including all mass dependencies has been completed in Refs. [31]. The 
genuine SUSY-QCD corrections at NLO are large, too, modifying the total gluon-fusion 
cross sections by up to ~ 100% depending on the MSSM parameters of the Higgs and 
stop/sbottom sectors. 

In this work we will present the soft-|-virtual-|-collinear gluon resummed results at 
N^LL for scalar Higgs boson production via gluon fusion including a consistent treatment 
of top and bottom mass effects with an extension to the inclusion of additional subleading 
collinear gluon effects. For the pseudoscalar Higgs boson we will present the analogous 
results up to the NNLL level, since the hxed-order N^LO result does not exist for this 
case. The paper is organized as follows. In Section 2 we will describe our method for the 
scalar and pseudoscalar Higgs bosons of the SM and MSSM. In Section 3 we will present 
results and discuss the comparison with previous calculations. In Section 4 we will hnally 
present our conclusions. 

2 SM Higgs Boson Production 

Higgs boson production via gluon fusion is mediated by top and to a lesser extent bottom 
and charm triangle loops at LO. Following the notation of Ref. [H] the production cross 
section of scalar Higgs boson production via gluon fusion gg ^ H can be cast into the 
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where the sum over i,j runs over all contributing initial state, i.e. only gluons at LO. The 
term denotes the corresponding parton-parton luminosity, Mh the Higgs boson mass, 
^R/F the renormalization/factorization scale, the strong coupling constant and Gf the 
Fermi constant. The variable tq is dehned as 
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where mq is the corresponding loop quark mass. The lower integration bound is given by 
tr = Ml/s. The integration kernel Gij can be computed as a perturbative expansion in 
the strong coupling constant. 
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where the LO and NLO expressions are given by 
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with the NLO quark mass dependence contained in the functions cr and {ij = 
99i9QiQQ)- Explicit results for these functions can be found in [20. In the limit of 

^Note that the results of [5] have to be divided by z due to our different normalization of the functions 
lj{z,TQ) {ij = gg,gq,qq). 
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heavy top quarks they approach the following simple expressions, 
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For the NLO results presented above we have used the notation (i = 0,1,...) 
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for the plus distributions and ^2 = 7r^/6. The Altarelli-Parisi splitting kernels are given 
by m 

Pggiz) = 

■1 
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The results of the NNLO pieces Gij of the coefficient function can be extracted from 
Refs. [miTB] and the Mellin transforms of their leading soft+virtual contributions are 
given explicitly in Ref. 


2.1 Soft and collinear Gluon Resummation 

In this work we will resum soft and collinear gluon effects up to all orders in the per¬ 
turbative expansion. This will be performed systematically in Mellin space. The Mellin 
moment of the gluon-fusion cross section is dehned as 
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The Mellin transformation can be inverted by means of the contour integral 
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where the value of the offset C has to be chosen such that all singularities of the N- 
moments are located on the left of the contour. Soft and collinear gluon effects arise in 
the limit z —)■ 1 that corresponds to the limit of large N in Mellin space. The leading 
contributions in Mellin space only appear in the gg initial state, while all other initial 
states are suppressed as 0{1/N) [18] as can be inferred explicitly from the NLO results 
presented in Eq. Following Refs. [T8K6] the leading contributions to the coefficient 
function Ggg in Mellin space are logarithmic in N and follow the perturbative expansion 
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where the second line indicates that the leading parts can be resummed in closed form. 
The all-order resummed expression can be cast into the form 
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where we included top mass effects up to the NLL level explicitly as will be discussed in 
following. The function Cgg contains all constant terms that originate from the 5(1 — z) 
terms of the explicit perturbative results and additional contributions emerging from the 
Mellin transformation. It develops a perturbative expansion in the strong coupling ag, 
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The leading logarithmic terms in N are contained in the exponential factor with the 
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where bo denotes the leading order term of the QCD beta function, 
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^Note that the soft and collinear limits of the mass-dependent functions dij{z,TQ) coincide with the 
heavy top limits presented in Eq. ([5]). 
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where Np is the number of active flavours that we choose as Np = 5 in the following, 
i.e. the top quark has been decoupled from the strong coupling and the PDFs. The 
term containing the function 

si;’(A) = |2^ + (1 - 2A) log(l - 2A)] (16) 

resums the leading logarithms in N, while the successive terms with {i > 2) cover 
the subleading, subsubleading, etc. contributions. The explicit expressions for and 
g^^ can be found in Refs. [T81IT7] and the one for g^^ in [H] where the four-loop cusp 
anomalous dimension has been approximated by employing Fade approximants. In the 
following we will use this approximate expression for g^^\ The NLO coefficient of the 
function Cgg is given by [IB] 

^99 = +6(7l + C2)-byclog^ 

hF 

= CHin) + 6 C 2 + log ^ (17) 

0 fip 

where we included the top quark mass dependence explicitly by using the function cuiTt) 
|20] . The NNLO coefficient (7® in the limit of heavy top quarks can be found in Ref. [T 8 ] . 
The N^LO term can be extracted from the explicit results given in [21] and can be found 
in Appendix E of [57]. In Ref. [12] the inclusion of additional terms of purely collinear 
origin in the resummation has been discussed and then implemented in the results of 
[I5|T8] by performing the replacement 

C™ ^ (18) 

This replacement reproduces all leading collinear logarithms of the form a”log^"^~^ N/N 
up to all orders. However, the subleading logarithms of this type are not covered by 
this replacement as can be inferred from the explicit NNLO expansion of the resummed 
expression in Ref. [18]. We have observed that the subleading logarithms of the type 
'-)E log^ N/N and log M]j/fXp log^ N/N can be reproduced by the extended replacement 



-)■ 


cli +6 


L 

N 


with the modihed logarithm 


L 


log 


Ne^^ fip 

Mh 


\ogN + 'yE - - log 


Ml 

/iF 


(19) 


( 20 ) 


This method can be extended to the next perturbative order. By performing the NNLO 
replacement 
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we reproduce the correct collinear logarithms log*iV/iV for i = 5,4 at N^LO, 

Gf^ = 36 log® N + 170.679 log® N + 744.849 log^ N + 1405.185 log^ N 
+ 2676.129 log^ N + 1897.141 log N + 1783.692 

+ ^{l08 log® iV +615.696 log^iV} + ... (22) 

where all numbers have been rounded at 0(10“^) and the additional contributions of 
the effective Lagrangian in the heavy top limit have not been taken into account as in 
Ref. [H]. These terms agree with the threshold expansion of the N^LO results of [21] and 
the logarithmic scale dependence terms up to the hrst two leading logarithms of 0{1/N). 
A further inclusion of subsubleading collinear logarithms requires a modihed systematic 
expansion of the Mellin transform of the resummed kernel in logA^ = log(iVe'^®) instead 
of log N and powers of N~^ which is beyond the scope of this work. It should be noted 
that our approach can be compared with the equivalent approach of Ref. [19] that tries to 
construct a second exponential resummation for the 0{1/N) logarithms to be added to 
the leading soft+virtual exponential. However, Ref. |19] does not apply their conjecture to 
the Higgs boson case so that a direct comparison of the methods is not possible presently. 

The numerical impact of the subleading collinear logarithms and top mass effects in 
the resummed expression ranges in the sub-per-mille range for the SM Higgs mass. The 
inclusion of top mass effects in the resummation can reach the per-cent level beyond NLO 
for large Higgs masses. An alternative implementation of resummation effects includes 
the full coefficient Cgg of Eq. (fT2|) in the resummed exponential with a careful expansion 
of the exponent up to all terms of 0{1/N) included in the analysis along the lines of 
Ref. [I2]. We have checked that this modihcation has an impact in the per-mille range 
on the hnal results for the inclusive SM Higgs cross sections. 


2.2 Mass Effects and Matching 


Using the resummed expression of the gluon-fusion cross section in Mellin space we sub¬ 
tract the corresponding hxed-order Mellin-space result up to NNLO in order to obtain 
the net effect of resummation beyond NNLO. This residual contribution has then been 
added to the NNLO result that has been obtained by including the full NNLO result in 
the limit of heavy top quarks for the pure top quark contributions and adding top mass 
effects and bottom (charm) contributions at NLO strictly. Furthermore we include only 
the top quark mass effects in the resummed cross section at the NLL level (convolved 
with NLO as and PDFs) and treat the bottom- and charm-induced parts at hxed NLO. 
Since the virtual coefficient of the bottom contributions behaves in the limit Mjj 3> 
as [9] {Ca = 3,Cf = 4/3) 
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if the bottom mass is renormalized on-shell, i.e. it contains large logarithms that are not 
resummed. The resummation of the Abelian part proportional to Cp has been performed 





in Ref. [50] up to the subleading logarithmic level. These logarithms are related to the 
Sudakov form factor at the virtual Hbb vertex that generates these large logarithmic 
contributions for far off-shell bottom quarks inside the corresponding loop contributions 
in the Abelian case. The resummation of the non-Abelian part proportional to the Casimir 
factor Ca has not been considered so far. This type of logarithmic contributions emerges 
from a different origin than the soft and collinear gluon effects discussed so far and is the 
main source of the very different size of QCD corrections to the bottom-loop contributions 
PllUmij . In order to obtain a reliable result for the bottom contributions a resummation 
of these types of logarithms is mandatory so that we do not include these contributions 
in our soft and collinear gluon resummation but treat them at hxed NLO. 

The complete cross section including soft/collinear gluon resummation effects can then 
be cast into the generic form 
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with (j/j denoting the LO cross section factor cxo of Eq. ([T]) including only the top quark 
contribution. The index \NNLO)^ in the second line indicates the fixed-order expansion 
of the resummed coefficient function G^gM in Mellin space up to NNLO while the index 
’’{NLOy denotes the perturbative expansion of the NLL resummed coefficient function 
&M^nll Mellin space up to NLO. The hrst integral has been convolved with N^LO 
and NNLO PDFs according to the discussion about the non-necessity of N^LO PDFs of 
Ref. [5l] and of resummed PDFs of Ref. [52] for the SM Higgs mass, while the second 
integral has been evaluated with NLO Og and PDFs consistently. The fixed-order NNLO 
cross section of the last term has been derived as 
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where denotes the full LO cross section including only top loops, the LO 

cross section including top and bottom/charm loops, (N)NLO K-factors 

obtained in the limit of heavy top quarks and the full NLO K-factor including 

top and bottom/charm loops. The NNLO parts have been derived with N^LO as and 
NNLO PDFs and the NLO terms with NLO as and PDFs consistently as implemented 
in the programs HIGLU [53] and SusHi [5l]. This implementation guarantees that top 
mass effects are treated at NLL level and bottom/charm contributions at fixed NLO 
respectively. 


2.3 Mellin Inversion 


As a final step we have to perform the Mellin inversion of the resummed and properly 
matched result in Mellin space. According to Eq. fl2T|) this requires the inclusion of 
the Mellin moments of the gluon densities /^(iV,in the integrand for the contour 
integration. This can be treated in different ways. Refs. [TSpH] have used a ht to the 
parton densities in Bjorken-x space for a fixed factorization scale and used this fit to 
determine the Mellin moments analytically in terms of a set of htted coefficients. We 
did not proceed along these lines. A second option will be to use the evolution program 
PEGASUS [55] in Mellin space that just requires the implementation of the input densities 
in terms of a predehned functional form. The latter does, however, not coincide with 
the input densities of all available global PDF fits. We have performed cross checks with 
PEGASUS for the MSTW08 PDFs [SB] at LO, NLO and NNLO(^. The method that we 
adopted in our numerical analysis is to implement the Mellin moments of the Bjorken-x 
PDFs by a numerical integration. Due to the fact that the 2 ;-integral of the hxed-order 
integral is extended to an infinite upper bound after integration [58] due to the presence of 
the Landau singularity at Ni = exp(l/[25oa<i(/^|;)]) that implies that the resummed kernel 
does not vanish for z > 1 but drops down very fast for larger values ol z. In order to 
increase the numerical stability of the Mellin inversion we have included four additional 
powers 1/(A^ — 1)^ in the resummed kernels so that the large N contributions for z > 1 
are sufficiently suppressed. By means of partial int^rations this can be translated to the 
convolution over derivatives of the gluon PDFs |59)u. 


a 


( 


2 ’ 2 


(27) 


with 


n{res) ( Ar. ( 2\ . ^h \ _ ^{res) ( at. ( 2\ . ^h \ / ( 1 \J _ 

5 ,.2 ’ ,.2 ) ^99 2 ’ ,.2 ) / ("^ 


2 ’ 2 


2 ’ 2 . 
h-ij h-F y 


^Differences between the original MSTW08 evolution and the results of PEGASUS at large Bjorken-x 
have been clarified after including an updated set of MSTW08 PDFs m- The impact on the gluon- 
fusion cross section is small, i.e. at the per-mille level both for small and large Higgs masses, and has 
been neglected in this analysis. 

^Note that the translation to the second derivative appearing in the Mellin-integral assumes that the 
gluon density and its first derivative vanish for x = 1. 
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The second derivative has been implemented by a equidistant three-point method 

fix) = /(^ + '■) - 2/W + /(^ - 1^) + (29, 

Using this method we can implement the original PDF hts in Bjorken-x space and obtain 
sufficient numerical stability if the contour is chosen according to the parametrization [59] 

N = C + (30) 


with the -l-(—) sign applied to the upper (lower) contour with respect to the real axis 
and the integration proceeds over the variable x. The off-set parameter has been chosen 
as C* = 2.5 and the angle 0 as Svr/d for z < 1 and vr/d for z > 1. This choice ensures 
that all relevant singularities in Mellin space are located to the left of the integration 
contour, but the Landau singularity to the right, i.e. the latter is excluded from the contour 
integral according to the minimal prescrition method of Ref. [58]. We have checked the 
independence of our results of variations of the two parameters C, 0 around the chosen 
values within the valid ranges, i.e. keeping C in the range 2 < U < Nl and the signs of 
the angle 0 and its hemispheres with respect to 7r/2. 


3 MSSM Higgs Boson Production 


In the MSSM all three neutral Higgs bosons are produced via gluon fusion gg —)■ h, H, A. 
In the following we will neglect the stop and sbottom loops in MSSM scalar Higgs boson 
production and will focus on the top and bottom-induced contributions that are modihed 
in the scalar case due to the additional mixing factors of the top and bottom Yukawa 
couplings compared to the SM Higgs case as given in Table [TJ This approximation works 
for heavy stops and sbottoms. Pseudoscalar Higgs boson production, however, requires a 
transformation of the SM Higgs results to the production of a CP-odd Higgs particle. The 
resummed exponential of Eq. flT^ is universal and thus unchanged in the pseudoscalar 
case (if Mh is replaced by the pseudoscalar Higgs mass Ma) so that only the coefficient 
Cgg is different from scalar Higgs boson production. The difference of this coefficient 
between the pseudoscalar and scalar Higgs cases is given bj@ 


gg 

gg 


c,. 

c. 


( 1 ) 

ggA 


( 2 ) 

ggA 

1939 

144 


- C’gUff = CAin) - Cnlt) 

^( 2 ) 

^gg,H 

9 >21,, /19 

+ STe + 6C2 - ~ 



Ml 

ml 


®Note that the NNLO expression differs from Ref. [55]. The (numerically minor) discrepancies have 
been clarified with the authors. 
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(31) 


(33-2Nf „ Ml 33-2 Nf, Ml 

where at NLO we again include the full top mass dependence in this coefficient along the 
lines of Ref. [20] . In the limit of heavy top quarks the difference of the virtual corrections 
at NLO approaches 

= C'L - ^ ^ (32) 

Since the gluon-fusion production cross section for pseudoscalar Higgs bosons is only 
known up to NNLO we reduce our resummation to the NNLL level, i.e. we do not include 
the function in the resummed exponent of Eq. ([HD. The inclusion of purely collinear 
logarithms proceeds along the same lines as for the SM Higgs case by implementing the 
replacements 


r'd) _v 

^IgA + (48 — iViJ’)— (33) 


with the extended logarithm 

~ , Ne^^giF , 1 , 

L = log - -= logiV + 7 e - - log 

Mh / Pi? 

These replacements resum the corresponding collinear logarithms at leading and sublead¬ 
ing level accordingly. 



4 Results 

We analyze Higgs boson production via gluon fusion at the LHC for a c.m. energy of 
y/s = 13 TeV for the SM Higgs boson and the neutral Higgs bosons of the MSSM at the 
N^LL level for the top-induced contributions of the scalar Higgs cases and at the NNLL 
level for the pseudoscalar case. As parton density functions (PDFs) we will primarily 
use the MSTW08 sets [56] with the strong coupling normalized to q;s(M|) = 0.12018 at 
NLO and to q;<j(M|) = 0.11707 (N)NNLO. The quark pole masses have been chosen as 
rrit = 172.5 GeV, = 4.75 GeV and rric = 1.40 GeV accordingly. For the determination 
of the PDF-t-CK^ uncertainties we have used the envelope method [60] with GTIO [HI] and 
NNPDF2.3 [62] PDF sets. For GTIO the strong coupling has been normalized consistently 
to as = 0.118 and for NNPDF2.3 to ag = 0.119, respectively. The uncertainty in the 
strong coupling constant has been adopted as ±0.002 around the corresponding central 
values of the used PDF^. 

®If other PDF sets as ABM12 [63] or HERAPDF1.5 [64] are included in this envelope the PDF-l-ag 
uncertainties will increase considerably with a major part originating from sizeable differences in the 
as fits at NNLO and different data sets included in the global fits. Moreover, the proper treatment of 
higher-twist effects in the global fits is an open aspect in this context. However, we have adopted the 
scheme used within the HXSWG, since an extended study of this particular issue is beyond the scope of 
our paper. 
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4.1 Standard Model 



Figure 1: Scale dependence of the SM Higgs production cross sections via gluon fusion for 
a Higgs mass Mh = 125 GeV and c.m. energy y/s = 13 TeV using MS TWOS PDFs. The 
renormalization and factorization scales are identified, i.e. fi = fiR = pip- 

For the SM Higgs boson we perform the numerical analysis for Higgs mass values 
around the measured Higgs mass of about 125 GeV. For the uncertainties we consider 
the scale dependence and the PDF+a* uncertainties. The scale dependence at LO, NLO, 
NNLO and N^LL is shown in Fig. [1] as a function of the identihed renormalization and 
factorization scales in units of the Higgs mass Mr- We observe a signihcant decrease of 
the scale dependence with increasing perturbative order signalizing a rather mild scale 
dependence at N^LL. However, this is caused by a signihcant cancellation between the 
renormalization and factorization scale dependences that develop opposite slopes [18] . 
For a consistent estimate of the theoretical uncertainties both scales have to be varied 
independently. The results of varying both scales by a factor of 2 up and down around 
our central scale choice /io = are shown in Table [2] without electroweak corrections. 

The numbers develop a much larger variation within the complete interval of indepen¬ 
dent variations so that identifying the renormalization and factorization scales leads to 
a signihcant underestimate of the related uncertainty. Taking the minimal and maximal 
values of Table [2] excluding the values for hf/tr > 2 and p.f/tr < | we derive a total 
scale dependence of [-1-4.4%, —5.3%], while for identihed scales we obtain the optimistic 
estimate of [4-4.4%,—0.1%]. Thus for a sophisticated determination of the theoretical 
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uncertainties due to the scale dependence both scales have to be varied independently 
contrary to what has been done in Ref. j22j . 


Hf/Mr 

1/4 

1/2 

1 

1/4 

46.08 pb 

42.92 pb 

39.35 pb 

1/2 

45.42 pb 

44.13 pb 

41.81 pb 

1 

44.98 pb 

45.81 pb 

44.94 pb 


Table 2: N^LL Higgs boson production cross sections via gluon fusion for different values 
of the renormalization and factorization scales Hr, without electroweak corrections for 
a Higgs mass Mr = 125 GeV and c.m. energy y/s = 13 TeV. 


Including electroweak corrections [2^ we present our predictions of the central cross 
sections for hve different Higgs mass values in Table |3] where we choose the MSTW08 PDFs 
[56] and add the corresponding PDF+a^ uncertainties. The total uncertainty is obtained 
by adding the scale and PDF+Os uncertainties linearly. In this way the total uncertainties 
amount to 0(10%). Additional uncertainties due to parametric uncertainties in the top 
and bottom/charm quark masses and missing quark mass effects beyond NLO are small, 
i.e. in the per-cent range and thus negligible. In Table H] we present our hnal predictions 


Mr [GeV] 

a{pp H + X) [pb] 

scale 

PDF-1-Q!s 

total 

124 

47.05 pb 

+4.5% 

-5.2% 

+3.7% 

-4.0% 

+8.2% 

-9.2% 

124.5 

46.72 pb 

+4.4% 

-5.3% 

+3.7% 

-4.0% 

+8.1% 

-9.3% 

125 

46.40 pb 

+4.4% 

-5.3% 

+3.7% 

-4.0% 

+8.1% 

-9.3% 

125.5 

46.06 pb 

+4.4% 

-5.2% 

+3.7% 

-4.0% 

+8.1% 

-9.2% 

126 

45.74 pb 

+4.4% 

-5.2% 

+3.7% 

-4.0% 

+8.1% 

-9.2% 


Table 3: N^LL Higgs boson production cross sections via gluon fusion for different values 
of the Higgs mass including the individual uncertainties due to the renormalization and 
factorization scale dependence and PDF+as uncertainties including electroweak correc¬ 
tions using MSTW08 PDFs for a c.m. energy yfs = 13 TeV. 


of the gluon fusion cross section for hve different Higgs masses with the central values 
and PDF-|-as obtained by the envelope method [60] involving MSTW08 [00|, CTIO [01] 
and NNPDF2.3 [02] PDFs. This results in a slight increase of the total uncertainties to 
the level of 0(15%). This situation will improve with the inclusion of more recent PDFs 
including LHC data in their global hts. 
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Mh [GeV] 

a{pp H X) [pb] 

scale 

PDF-l-a^ 

total 

124 

48.37 pb 

+4.5% 

-5.2% 

±8.1% 

+12.6% 

-13.3% 

124.5 

48.00 pb 

+4.4% 

-5.3% 

±8.1% 

+12.5% 

-13.4% 

125 

47.63 pb 

+4.4% 

-5.3% 

±8.1% 

+12.5% 

-13.4% 

125.5 

47.28 pb 

+4.4% 

-5.2% 

±8.2% 

+12.6% 

-13.4% 

126 

46.94 pb 

+4.4% 

-5.2% 

±8.3% 

+12.7% 

-13.5% 


Table 4: N^LL Higgs boson production cross sections via gluon fusion for different values 
of the Higgs mass including the individual and total uncertainties due to the renormaliza¬ 
tion and factorization scale dependence and PDF-hag uncertainties including electroweak 
corrections using the envelope of MS TWOS CTIO l6lf and NNPDF2.3 J6^ PDFs 

for a c.m. energy s/s = 13 TeV. 


Finally, we compare our results with the N^LO predictions of Ref. [22] in Table |5l This 
table shows the numbers of Ref. [22| in the hrst column which have been obtained in the 
limit of heavy top quarks for the QCD corrected cross sections at all perturbative orders 
and neglecting the bottom and charm loops. The renormalization and factorization scales 
have been identihed for the derivation of the scale dependence given in per-cent attached to 
each number. The second column labelled ’SVCqo’ depicts our prediction in the heavy top 
limit by using the approximate soft-|-virtual-|-colhnear expansion of our resummed kernel 
for the N^LO piece added to the full NNLO result in the same limit. These approximate 
hxed-order results agree with the explicit numbers of Ref. [22] within 0.2% for /i/j = 
fip = Mh/2 and within 1.5% for fiR = fip = Mh and develop a similar scale dependence. 
Our corresponding result at N^LL is shown in the third column labelled ’N^LLqo’ in 
the limit of heavy top quarks and omitting the bottom and charm contributions. Our 
scale dependence for identihed scales is of similar size as the one of the N^LO results. 
However, this scale dependence does not constitute a reliable estimate of the theoretical 
uncertainties as discussed before. The comparison of the second and third column shows 
the effect of resummation beyond N^LO. It is clearly visible that the resummation effects 
range at the per-mille level for the scale choice p-R = p-F = Mh/2, while resummation 
provides a 2% contribution beyond N^LO for the scale choice Pr = Pf = Mh- The 
latter effect is of the order of the uncertainties related to the scale dependence. The 
fourth column presents our full N^LL results including top mass effects up to NLL and 
bottom/charm loops up to NLO. The value for our central scale choice Pr = Pf = Mh/2, 
agrees with the number of Ref. [22] within 0.3% by accident. The last column shows 
our numbers including the electroweak corrections of Ref. [27]. Our hnal results deviate 
by a few per cent from the numbers presented in [241125] . since the hrst analysis works 
at approximate hxed N^LO and does neither include top mass effects beyond LO nor 
bottom/charm loops and in both analyses electroweak corrections have not been taken 


15 















p/Mh 

Ref. [22] 

SVGoo 

N^LLo, 

massive N^LL 

with elw. corr. 

1/2 

44.31«;g; pb 

44.22+t^l pb 

44.15l^;g pb 

44,13«;« pb 

46.40l^:g pb 

1 

ii.lYlil pb 

43.77^;S pb 

44.72^5;^ pb 

44.94+®;®^“ pb 

47.251?;^! pb 


Table 5: Higgs boson production cross sections via gluon fusion for different values of the 
identified renormalization and factorization scales p. = fiR = fvp with and without guark 
mass effects and electroweak corrections in comparison to the N^LO results of Ref. J2^ 
using MSTW08 PDFs for a c.m. energy ^/s = 13 TeV. The percentage errors provide the 
relative scale dependence for identified renormalization and factorization scales p = Pr = 
Pf varied by a factor of 2 up and down around the central scale choice. The second and 
third columns of cross sections in the heavy top limit at all perturbative orders, labelled 
’SVCoo’ and ’N^LLoo’, exclude bottom/charm contributions. 


into account. Moreover, small differences can arise due to different implementations of 
the PDFs in the resummation framework in comparison with the second analysis [25] . 

4.2 Minimal supersymmetric extension 

For neutral MSSM Higgs boson production via gluon fusion we have adopted the 
scenario of Ref. [65] which is dehned as 

Msusy = 1 TeV, p = M 2 = 200 GeV, Xt = I.QMsfsyi 'oOg = 1.5 TeV, A^, = Aj. = A^ (35) 

where Msusy denotes the SUSY-breaking sfermion mass scale of the third generation, 
p the higgsino mass parameter, M 2 the wino mass parameter, Xt = At — ptgfi the stop 
mixing parameter, the gluino mass and At^^r the trilinear SUSY-breaking couplings 
of the sfermions and Higgs helds. We use the RG-improved two-loop expressions for the 
Higgs masses and couplings of Ref. [66] which yield predictions for the Higgs boson masses 
that agree with the diagrammatic calculations of Ref. [30] within 3-4% in general. Thus 
the leading one- and two-loop corrections have been included in the Higgs masses and the 
effective mixing angle a. 

In Fig. [2] we show the scalar and pseudoscalar production cross sections as functions of 
the corresponding Higgs masses at LO, NLO, NNLO and (N)NNLL for the pseudoscalar 
(scalar) MSSM Higgs bosons. Squark loops and genuine SUSY-QGD corrections have 
been neglected in this work. They can be added in a form factorized from the resummation 
effects which is left for future work. The resummation effects beyond NNLO amount to 
about 5% for tgfi = 3, where the top loop contributions are dominant, while for large 
values of tgfi = 30 the effects are small, since the bottom loops dominate so that the 
accuracy of the cross section is of hxed NLO according to our setup as discussed in 
Section [2?2l This can be inferred more clearly from Fig. [3] which shows the corresponding 
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K-factors for scalar and pseudoscalar Higgs boson production via gluon fusion. The K- 
factors are defined as the ratios of the NLO, NNLO and (N)NNLL cross sections to the 
LO prediction where each order is conistently evaluated with the corresponding PDF and 
as choices. For large values of tg/3 = 30 it is obvious from this figure that there is only a 
tiny effect beyond NLO as expected due to the dominance of the bottom loops that are 
purely treated at hxed NLO. The bumps and spikes for Mh/a ~ 2mt are related to the ti 
threshold that generates a Coulomb singularity for the pseudoscalar case already at NLO 
PI3^ . The latter is regularized by taking into account the hnite width of the virtual top 
quarks which, however, is beyond the scope of our work. The results indicate that the 
QCD uncertainties reduce to the level of ~ 10% after including the resummation effects 
in regions of top-loop dominance. The dominant uncertainties in these regions within the 
MSSM will arise from genuine SUSY effects. 

5 Conclusions 

In this work we have performed an analysis of soft and collinear gluon resummation effects 
in SM and MSSM Higgs boson production via gluon fusion at the LHC. For scalar Higgs 
boson production the results have been analyzed at the N^LL level, while pseudoscalar 
Higgs boson production can only be extended to the NNLL level with the present state- 
of-the-art calculations. We have implemented top mass effects up to the NLL level by 
using the full mass dependence of the hnite part of the virtual corrections according to 
PD]. The bottom and charm loops, however, have been treated at hxed NLO rigorously 
since there are additional non-resummed double logarithmic contributions that have a 
large impact on the size of their contribution. We have discussed an alternative extension 
of the formerly used approach for the inclusion of collinear gluon ehects at 0{1/N) in 
Mellin space to the subleading logarithmic level as a conjectur^. We have been able 
to reproduce the leading and subleading N^LO terms at 0{\/N) which is a non-trivial 
cross check of our method. The impact of resummation ehects on the total cross sections 
reaches a size of a few per cent beyond the hxed-order calculations, while the ehect of 
including top mass ehects and subleading collinear logarithms in the resummation ranges 
at the per-mille level or below. 

Acknowledgements. We are indebted to C. Anastasiou, J. Bliimlein, D. de Florian, 
M. Grazzini, S. Moch, J. Rojo, R. Thorne, A. Vogt and J. Zurita for useful discussions 
and helpful clarihcations. We are grateful to R. Thorne for providing us with an updated 
version of the MSTW08 PDFs. Ths work is supported in part by the Research Executive 
Agency (REA) of the European Union under the Grant Agreement PITN-GA-2012-316704 
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^The recently developed next-to-eikonal approach [67] may provide a basis for a factorization proof of 
logarithmic 0{\/N) terms. 
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[GeV] 

Figure 2: MSSM Higgs production cross sections via gluon fusion at the LHC for a c.m. en¬ 
ergy y/s = 13 TeV and two values of tgfl = 3,30 using MSTW08 PDFs. For the MSSM 
the scanario fSBf has been adopted and sguark loops as well as genuine SUSY-QCD 

corrections have been neglected. 
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Figure 3: K-factors of the MSSM Higgs production cross sections via gluon fusion at the 
LHC for a c.m. energy y/s = 13 TeV and two values of tg(3 = 3,30 using MSTW08 
PDFs. For the MSSM the scenario has been adopted and sguark loops as well 

as genuine SUSY-QCD corrections have been neglected. 
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